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Hole maximum density droplets of an antidot in strong magnetic fields 
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We investigate a quantum antidot in the integer quantum Hall regime (the filling factor is two) 
by using a Hartree-Fock approach and by transforming the electron antidot into a system which 
confines holes via an electron-hole transformation. We find that its ground state is the maximum 
density droplet of holes in certain parameter ranges. The competition between electron-electron 
interactions and the confinement potential governs the properties of the hole droplet such as its spin 
configuration. The ground-state transitions between the droplets with different spin configurations 
occur as magnetic field varies. For a bell-shape antidot containing about 300 holes, the features 
of the transitions are in good agreement with the predictions of a recently proposed capacitive 
interaction model for antidots as well as recent experimental observations. We show this agreement 
by obtaining the parameters of the capacitive interaction model from the Hartree-Fock results. An 
inverse parabolic antidot is also studied. Its ground-state transitions, however, display different 
magnetic-field dependence from that of a bell-shape antidot. Our study demonstrates that the 
shape of antidot potential affects its physical properties significantly. 

PACS numbers: 72.15.Qm, 73.23.Hk, 73.43.-f 
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I. INTRODUCTION 



A quantum antidot has been extensively investigated 



.ntidot has been extensively mvcstiga 
experimentally [1 liaji iijl |, ft M El lU and 

theoretically [1 11 111 111 El 13 last decade. It is a 
potential hill in two dimensional electron gas (2DEG) 
systems. In zero magnetic field, it is a simple repulsive 
potential and acts as a scattering center for electrons. 
In this sense, it is opposite to a quantum dot which 
confines electrons. When a strong magnetic field is ap- 
plied perpendicular to 2DEG, the antidot has its own 
electronic "edge" structures, which correspond to clas- 
sical skipping orbits around the antidot resulting from 




FIG. 1: A quantum Hall antidot with localized antidot states 
and extended edge channels. The localized states are weakly 
coupled to the extended channels, as indicated by dashed and 
dotted line. The small arrows show the spin direction of the 
states. The local filling factor around the antidot is two. 
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the Lorentz force. These localized antidot structures 
can be experimentally studied by measuring conductance 
when they are weakly coupled to extended edge channels 
propagating along the boundary of 2DEG. In the integer 
quantum Hall regime, the measured conductance exhibits 
interesting Aharonov-Bohm oscillations H 0, II U ^3 > 
which can not be understood within a single-particle pic- 
ture. For exampleSjthe oscillations are accompanied by 
the charging effect la, SB, nontrivial h/(2e) Aharonov- 
Bohm oscillations HE3l, and/or Kondo-like signatures 
E3- These observations indicate that electron-electron 
interactions can play an important role in the antidot 
system. 

However, there have been few theoretical works on the 
interaction effects in the antidot system. Very recently, a 
phenomcnological capacitive interaction model has been 
proposed EH to explain the experimental results. This 
model is based on the capacitive couplings between local- 
ized excess charges, which are formed around the antidot 
due to magnetic flux quantization. The capacitive inter- 
action of the excess charges results in Coulomb block- 
ade and tunnelings are allowed only under certain con- 
ditions. The main result of the above mentioned work 
was that the usual resonant tunnelings are accompa- 
nied by Kondo resonances, hence leading to nontrivial 
Aharonov-Bohm oscillations. These predictions are in 
qualitatively good agreement with the experimental ob- 
servations H E3- Based on a Hartree-Fock approach 
and a particle-hole transformation, it was also suggested 
|l5lll6j that holes inside an antidot can form a maximum 
density droplet (MDD) El El El HI in the ground 
state within some parameter ranges and that the tran- 
sitions between MDDs may lead to Kondo effects, sup- 
porting the capacitive interaction model. However, the 
tested antidot was so small (it has about 50 holes) that 



2 



the transitions did not occur periodically with varying 
magnetic field, in contrast to the experimental data. 

In this paper, we develop a microscopic Hartree-Fock 
approach and apply it to a large-size antidot containing 
about 300 holes. Our approach is based on an electron- 
hole transformation, where an antidot potential of elec- 
trons is transformed to a confinement potential of holes. 
As in the experiment |9|, |l2|, we consider the antidot 
states formed by electrons with spin up and down in the 
lowest Landau level (i.e., the local filling factor is two 
around the antidot as in Fig. We test two kinds of 
antidot potentials: bell-shape and inverse parabolic. 

For both potentials, the antidot ground states are 
found to be MDDs of holes in certain parameter ranges 
(see Fig0). For a given magnetic field, the spin config- 
uration of MDD ground states (the size and the spatial 
splitting between spin-up and -down edges of the droplet) 
is determined by the competition |l7j| between electron- 
electron interactions and the confinement potential of 
holes: droplets with larger size are favored by weaker 
confinement and stronger electron-electron interactions. 
As magnetic field varies, the relative magnitude of these 
two competing factors changes so that the transitions of 
MDD ground states can take place. In general, there 
can be three types of transition: spin-up, spin-down, and 
spin-flip transitions, which manifest themselves via spin- 
up electron normal resonance, spin-down electron normal 
resonance, and Kondo resonance, respectively. 

For a bell-shape antidot potential, we find spin-down 
(Fig. [2J) and spin-flip transitions (Fig. |3J) in some ranges 
of magnetic field, while we do not find spin-up transi- 
tions. In the spin-flip transitions, the number of spin- 
down (spin-up) holes decreases (increases) by one as mag- 
netic field becomes stronger. A series of these transitions 
is obtained as a function of magnetic field in Fig. The 
features of the transitions are in agreement with the pre- 
dictions of the capacitive interaction model |15| , and thus 
one can explain the experimental observation [T^ | from 
them. From the Hartree-Fock result of the transitions, we 
obtain the parameters of the capacitive interaction model 
and study the variation of excess charges as a function of 
magnetic field. 

For an inverse parabolic antidot potential, we find 
all the three types of transitions in some parameter 
range. The properties of the spin-flip transition are dif- 
ferent from those of the bell-shape potential; in this case, 
the number of spin-down (spin-up) holes increases (de- 
creases) by one as magnetic field increases. Moreover, the 
spin-flip transitions are found to appear more frequently 
than the spin-down and spin-up transitions. As a result, 
the spatial splitting between spin-up and -down edges 
of MDDs becomes larger for stronger magnetic fields. 
This behavior differs from the experimental situation of 
Kataoka et ai, where the splitting is expected to be a pe- 
riodic function of magnetic field. It would be interesting 
to investigate inverse parabolic antidots experimentally 
and compare the obtained results with our theoretical 
predictions. 
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FIG. 2: Schematic diagram of particle densities, as a func- 
tion of distance from antidot center, in a hole MDD where Nf 
spin-up holes occupy single-particle states with angular mo- 
mentum m = 0, 1, • • • , iVf — 1 and Ni spin-down holes with 
m = 0, 1, • • • , Ni — 1. Single-particle hole states with spin a 
are empty for m > N a — 1. Note that a single-particle state 
with smaller m is located at smaller distance from the cen- 
ter. The changes in electron and hole densities around the 
antidot are also shown when a spin-down MDD transition 
\Nf,Ni) — ► \Nf,Ni + 1) occurs. In this process, a spin-down 
electron tunnels out of the antidot (indicated as thin arrows) 
and thus the total hole spin decreases. 
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FIG. 3: Same diagram as in Fig. H but for the spin-flip 
transition \Nf, iVj,) — > |iVj + 1, JVj_ — 1). In this process, a co- 
tunneling event (see thin arrows) takes place, where a spin-up 
electron moves out of the antidot while a spin-down electron 
moves in, and thus the total electron (hole) spin decreases 
(increases). 
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Hamiltonian: 
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B[T] 
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FIG. 4: Chemical potential, /j,n = -Ejv+i — En, vs magnetic 
field B for a bell-shape antidot potential. Different values of 
TV = iV T + 7V| G [368, 384] are used. The horizontal dotted 
line represents hole Fermi energy, which can be rather dif- 
ferent from electron Fermi energy. Diamonds represent the 
spin-flip transition \N^,Ni) -» |/V t + 1,/V; - 1) as B in- 
creases, while vertical jumps show the spin-down transition 
| AT-j- , ATj, ) — » \Nt;,Ni + 1). Following the zigzag solid line, 
the MDD ground state evolves as \N r ,N L ) = 1 175, 193) -> 
|175, 194) -» |175,195> 1 176, 194) -> 176, 195) -> 

1 176, 196) [177,195) -> [177,196) -> 177, 197) => 

1 178, 196) -» 1 178, 197) -> [178, 198) with increasing B. Here 
— > and => indicate, respectively, the spin-down and spin-flip 
transitions. The parameters of this antidot can be found in 
Sec. VB. 



This paper is organized as follows. In Sec. II our model 
Hamiltonian for the antidot is given. This Hamiltonian 
is changed into a hole Hamiltonian via an electron-hole 
transformation in Sec. III. Within a Hartree-Fock ap- 
proach, we study the stability of MDDs in Sec. IV and 
the MDDs of a bell-shape antidot in Sec. V. In Sec. VI 
the properties of the MDDs of the bell-shape antidot are 
shown to be in good agreement with the prediction of the 
capacitive interaction model. In Sec. VII we investigate 
the MDDs of an inverse parabolic antidot. Concluding 
remarks are given in Sec. VIII. 



II. MODEL HAMILTONIAN OF ANTIDOT 

We consider a 2DEG around an antidot in the pres- 
ence of a strong perpendicular magnetic field B along 
the z-axis. Following experiments 0, 0], the local 
filling factor around the antidot is chosen to be two. 
Thus, the antidot states can be assumed to be formed 
by spin-up and -down electrons in the lowest Landau 
level. In the symmetric gauge, the single electron wave 
function (f> m (r) is labeled by the quantum number m, 
the z-component of the angular momentu m. More e x- 
plicitly, (f> m (r) = (z/e) m exp(-\z\ 2 /4e 2 )/{\/2 m + 1 Trm\e), 
where z — x + iy is the complex coordinate of the 2D 
plane and 1(B) = yJhc/eB = 2.56 x 1CT 6 / ^/B[T][cra] 
is the magnetic length. We can write the antidot model 



H = EE^h 
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(i) 



where W(m) is the antidot potential energy [see Eqs. 10 
|3J| below], V is the Coulomb interaction, and c^ ma cre- 
ates an electron in the state <j) m (f) with spin a. The 
term of V™ n comes from the neutralizing positive back- 
ground charge around the antidot. Since there are as 
many ions as the total number of electrons, both spin-up 
and -down, we have V™ n = 2^™° (mm'\V\mm') . By 
definition V™ n > 0. The angular momentum conserva- 
tion yields mi' + 7712' = mi + m^- In our numerical work 
we use electron states up to a cutoff value m c , which is 
chosen sufficiently large. 

The following Coulomb matrix elements |2l| are used: 



(m + p, n\V\m, n - 
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This expression allows us to perform Hartree-Fock calcu- 
lations even when the total number of holes inside antidot 
is quite large (more than 300). 

We consider two types of antidot potential in this pa- 
per. The first type is an inverse parabolic potential 



W(r) 



\hw c 

constant, 



\m*VL 2 r 2 , r < r s 



r > r s 



(2) 



Beyond r > r s the potential is flat as a function of r. 
The second type is a bell-shape potential 



W(r) = { B 

constant, 



c 

77 , 



r <r t 

r t < r < r s 

r > r s . 



(3) 



Here tu c = \e\B/(m*c), electron charge e < 0, and 
m* = 0.067m e for GaAs. In the interval r < r t the 
potential is inverse parabolic, while in the next interval 
rt < r < r s the curvature changes sign. The matrix ele- 
ments W(m) = (m\W(r)\m) can be approximately writ- 
ten as follows when r t 3> I and r s — r t S> t: 



rT . 1 IliWc-fi fiVfm + 1), m < m s 

W(m) = < 2 y ' 

1 constant, 



m > m, 



(4) 
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for the inverse parabolic potential, while 



W(m) 




m*Cl 2 £ (m + 1), m < m t 



constant, 



m t < m < m s (5) 
m > m s 



for the bell-shape potential. Here mi — r 2 /(2l 2 ) — 1, 
where i = s, t. In order to make sure that the bell- 
shape potential is continuous at m = m t , we set Wo = 
fuj c /2 - m*n 2 £ 2 (m t + 1) - C/(2l 2 m t ). Note that W(m) 
is monotonously decreasing with m. 

We remark that the values of r t and r s are fixed by 
the shape of W(r), while mt and m s are magnetic-field 
dependent. But we can use constant m t and m s since, 
for the narrow range of magnetic field (SB « 0.02 T) of 
interest in this work, the change Srrii ~ miSB / (2B) is of 
the order of one, much smaller than m t and m s . In our 
bell-shape antidot, both positions of spin- up and -down 
edges of MDDs are in the interval of r t < r < r s so that 
they feel 1/r 2 potential. Note that m s is different from 
m c , which is the cutoff for the single particle states. In 
our numerical work we choose m c = 400 and m s = 300. 
The typical value of mt is between 100 and 200. 



III. ELECTRON-HOLE TRANSFORMATION 

An electron antidot system is an open geometry prob- 
lem and often requires heavy calculations to compute 
its physical properties. Such a difficulty can be avoided 
by transforming an electron antidot system to a system 
which confines holes since the transformed system con- 
tains only a finite number of holes. Such transformation 
is described in this section. 

We consider a particle-hole transformation 0, [2^J| of 
the type c m<a — > h^ m a and c] n (J — > h m , a . The term 
representing the interaction with the positive background 
in Eq. Jpl transforms into 

m c m c m c 

/ , / t v m ^ma^ma — & / J / J v m ll ma ll ma -± / j v m , 
am am m 

where = (mm'\V\mm') . The total Hamiltonian 
in Eq. JQl can be rewritten as 

H = ^ £ m h mt h m <f + e m h mi h m x 

m in 

+ 2 SflB J2 ( h U km1 i ~ k lnl h ™l 
m 

^ m c 

+ 2 ^ ^ (m^m'^Vlm^) 

a\a2 mi'm^'m^m^ 

X h , h 

A "m 2 i72"miiTi"nilVi"m! , ii! 



where V* — Ylm° (mm'\V\m'm). The effective single 
hole energy is 



-W(m) + V r 



x 



(7) 



The first term of Eq. J2J is the confinement energy com- 
ing from W(m) after the transformation, while the sec- 
ond term represents the change in exchange energy when 
an electron with angular momentum m disappears. A 
Hartree term is absent in e m , since it is canceled by 
the interaction between a hole with the quantum num- 
ber m and the positive background charge. So in this 
transformed Hamiltonian of hole systems the background 
charge term is absent. Note that the Zeeman term in Eq. 
© has the opposite sign to the corresponding term in 
Eq. P. 

There is a simple check of this result. If there are 
zero holes, according to Eq. ©, the total energy is 
Et = 2 Em W(m) - 2 YTm Vm ~ Em V m ■ This energy 
should be equal to the total energy of an electron antidot 
system in which electrons are occupied from to m c ; in 
this case the Hartree-Fock theory exactly gives the total 
energy, which has four contributions of the confinement 
energy (2^™ c W{m)), interaction energy with the pos- 
itive background (-4£™ c Vp, Hartree (2£™ c ^), 
and exchange energies (— X]™ c V*). The constant terms 
(Et) in the hamiltonian of Eq. © will be ignored here- 
after. 

When the cutoff value m c is larger than the total num- 
ber of holes, it is a good approximation to treat V* in 
Eq. J7J) as a constant Then, e m is a monotonously 
increasing function of m. For the bell-shape potential, 
for example, de m /dm ~ —dW(m)/dm = C / (2m 2 £ 2 ) > 
in the interval mt < m < m s [see Eq. ©]. 

Note that the magnetic field dependence of the hole 
confinement potential is 



W(m) 



1/B, for parabolic 

B, for bell-shape (~ 1/r 2 region). 



(8) 



X 
m ■ 



(6) 



On the other hand, the interaction energy scale e 2 /(e£) 
is proportional to \f~B. Therefore, the hole con- 
finement potential relative to the interaction energy, 
—W(m)/[e 2 /(e£)], decreases with increasing B for an in- 
verse parabolic antidot, while it increases in the 1/r 2 
potential region of a bell-shape antidot. 



IV. HOLE MAXIMUM DENSITY DROPLETS 

In electron quantum dots, MDDs are exact ground 
states when a strong magnetic field is applied and the 
confinement potential is strong enough. The reason is as 
follows: A MDD is an eigenstate of L z , the ^-component 
of the total angular momentum. In fact, it has the small- 
est possible eigenvalue of L z for a given number of elec- 
trons and there are no other states in the Hilbert space 
with the same eigenvalue. If one chooses a rotation- 
ally symmetric potential, a MDD is also an eigenstate 
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FIG. 5: The boundary (critical magnetic fields) of sta- 
ble and unstable MDDs of a bell-shape antidot. The fol- 
lowing parameters are used: HQ = 1.5 meV, C/(2£ 2 ) — 
2100^/B (Tesla) x e 2 /(<0, and m t = 118. 



of the Hamiltonian. Since the mean radius of the single 
electron wave packet increases with the ^-component of 
the angular momentum, the confinement potential favors 
small values of L z . Therefore, a MDD is certainly the ex- 
act ground state if the confinement potential is infinitely 
strong. It must also remain so in a certain parameter 
range as long as the potential is strong enough. The 
properties of a MDD were investigated by exact diag- 
onalization for a small dot ,17]. For a dot with about 
50 electrons, the properties of spin-polarized MDDs and 
their instability were also studied 0, 0] using Hartree- 
Fock approach and exact diagonalization. Experimental 
investigation of electron MDDs were reported by several 
groups [2fT |. 

Our antidot problem becomes similar to the quantum 
dot case after the particle-hole transformation described 
in the last section. Like electron MDDs, hole MDDs have 
a single-Slater-determinant form, 



I^T^i) = <- U ' ' ' A_i,t ' • ^otI°>- ( 9 ) 



The total number of holes is N = + N^, and N± is 
equal to or larger than iVf due to the Zeeman energy. 

The hole MDD is excellently described by the Hartree- 
Fock approach. Its total Hartree-Fock energy is 
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FIG. 6: The renormalized Hartree-Fock single hole energies 
£„ct °f a MDD are plotted for both spin up (solid) and down 
(dashed line). This MDD is stable as < £jv ct -ict f° r 

all m < JV CT — 1 and for both spins. We choose parameters 
Ml = 1.5 meV, C/(2£ 2 ) = 2100,/B (Tesla) x e 2 /(0, m t = 
118, N = 355, and B = 2.13T. The last occupied states 
(m = Ncr — 1) are marked with vertical lines. Note that 
iV T = 152 and N l = 203. 



Zeeman energies, respectively. Here 

JV T -1 Nf-l N l -1 N I 
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(-W(m) + V X )+Y, (-W(m) + V x ) 



In the above expressions, we have used the definitions 
Vmm' = {mTn'\V\mm') and V x m , = (mm'\V\m'm). 

From Eq. I|1U|) . one can define the renormalized single 
hole energy, which includes the Hartree and exchange 
self-energy corrections 



-HF 



W(m) + V„ X +J2 V mm' 
m' 

N--1 N a -1 
E V "m'~ E V X m , + gl lB S<T , (11) 



E HF (N h N l 



E 



H 



E 



X 



■ Ez + E, 



c, 



(10) 



where Eh, Ex, and Ez are the Hartree, exchange, and 



where Sf = 1/2 and sj, = —1/2. This renormalized single 
hole energy is useful for studying the stability of MDD 
states. A MDD state will be a stable ground state if 
the occupied single hole states satisfy s^a < £ (a£_i) ct 
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FIG. 7: Same curves as in Fig. ©but at B = 2.205T. This 



MDD is unstable since e 
that iV T =153 and N t = 202 



for m f» 125. Note 



for all m < A^ — 1 and for both spins. For a given N, 
this condition is satisfied, i.e., MDD states are stable, 
below some critical magnetic fields. A phase boundary 
between stable and unstable states is displayed in Fig. 
[3] In calculating the phase boundary for a given value of 
N, we calculate the ground configuration |iV-f,JVj,) and 
£ ma with varying B, and determine the critical value 
of B, where the above mentioned stability condition is 
violated. The examples of stable and unstable MDDs 
are plotted in Figs. [Uandd respectively. 

The physics determining the spin configuration 
(ATf, Ni) of a MDD ground state is as follows: There is a 
competition between the Coulomb energy (Eh + Ex) 
and the confinement energy (Eq)- If the confinement en- 
ergy is strong, the total energy is minimized by making 
the confinement energy small, i.e., by making the droplet 
size small. As a result, for a given N, the configuration 
with smaller N± — is favored by stronger confinement, 
since the droplet size is determined by N± . On the other 
hand, if the confinement energy is weak, the total en- 
ergy can be minimized by making the Coulomb energy 
smaller, i.e., by making the droplet size larger. Thus the 
degree of spatial splitting between spin-up and -down 
edges (or N± — ATf) of MDDs depends on the relative 
strength of the confinement energy and electron-electron 
interaction. 



HARTREE-FOCK RESULTS OF 
BELL-SHAPE ANTIDOTS 



In this section, we study the bell-shape antidot with 
the potential of Eq. J3J. Before discussing the Hartree- 
Fock result of its hole ground states, it is instructive to 
consider possible transitions of the hole ground states. 
These transitions can occur since the competition be- 
tween the Coulomb energy and the confinement energy 
of MDDs varies with magnetic field. They correspond 



to resonant tunneling processes 0] of the antidot when 
the antidot states are weakly coupled to extended edge 
channels (see Fig. QJ. 



A. possible ground state transitions 

In general, there exist three kinds of transitions be- 
tween MDD ground states, which are spin-up, spin-down, 
and spin- flip transitions. Within the Hartree-Fock ap- 



proximation, the spin-up transition of \N+ 
l.N 



N 



|JV T ± 



I) occurs at the magnetic fields where the degeneracy 
E HF (N h N l )=E HF (N 7 ±l,N l ) (12) 



is satisfied. The plus and minus signs refer to tunneling in 
and out of a hole, respectively. Similarly, the spin-down 
transition occurs when 



Ehf(N^, 2Vj_) = Ehf{N], Ni ± 1). 



(13) 



In electron language, the spin-up (-down) transitions cor- 
respond to normal resonant tunnelings of spin-up (-down) 
electrons fl5| . The normal spin-down resonant tunnel- 
ing, for example, is illustrated in Fig. [3 With increasing 
magnetic field, the direction of these transitions is cither 
\N h Ni) -> |AT T + 1,AT;) or |AT T ,AT ; ) \N h N t + 1) so 
that the total number (N = Nf + JVj_) of holes becomes 
larger. This increase of the total hole number can be un- 
derstood from the fact that the effective hole potential 
— W(m) decreases with magnetic field while the Fermi 
energy does not change. On the other hand, the spin-flip 
transitions occur at the magnetic field where the degen- 
eracy 

E HF (N h N l )=E HF (NT±l,N l Tl) (14) 

is satisfied. These processes can cause the Kondo 
resonance [2j|, [24| in the antidot system when both the 
spin-up and -down parts of MDDs are weakly coupled to 
the corresponding extended edge channels; this resonance 
is accompanied by the cotunneling processes where an 
electron (hole) tunnels into antidot and another electron 
(hole) with the opposite spin tunnels out via a virtual 
state |N T ± 1, Ni) or \N h N t t 1). 

We note that as magnetic field increases, the spin-flip 
transitions occur as either |JV-f,iVj) — > | JV-] — hi, iVj_ — 1) (see 
Fig. or |jVf, Ni) -> |A^ T - 1, JV; + 1>, depending on the 
detailed shape of antidot potential such as its curvature 
near MDD edges. For the bell-shape antidot (see Sec. 
VB), we find the former direction of transition, while the 
latter for the inverse parabolic antidot (Sec. VII). 



B. Hartree-Fock results 

In this subsection, we perform Hartree-Fock calcula- 
tions for the bell-shape antidot with more than 300 holes 
and study the transition of MDD ground states. We find 
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B [T] 

FIG. 8: Same curves as in Fig. [I]but for a weaker bell-shape 
potential. Twenty one different values of JV £ [350, 370] are 
used. Following the topmost zigzag solid line, for example, 
the MDD ground state evolves as |JVj,JVj) = 1 128, 233) -> 
|128, 234) |129, 233) -» |129, 234) -> |129, 235) => 

|130, 234) -» |130, 235) -» |130, 236) => |131, 235) -> 
|131, 236) -> |131, 237) |132, 236) with increasing B. 
Similarly, in the fourth zigzag solid line, |JVj,JVj) evolves 
as 1 127, 231) -» 127, 232) => |128, 231) -v |128, 232) -> 
1 128, 233) [129,232) -> [129,233) -> 129, 234) => 

1 130, 233). Note that both of 7V T and JVj increase with B. 

the spin-flip and spin-down transitions in certain param- 
eter ranges. 

From the Hamiltonian of Eq. 1)11)0. we find the energy 
of MDD ground state |JVf,JVj) with varying magnetic 
field B for a given total number of holes JV = iV-j- + 
JVj. In Fig. [SI the chemical potential /ijv = -Ejv+i — 
-EjVi which is the energy difference of ground states, is 
plotted as a function of B. In thi s calculatio n, we use 
Ml = 1.5 meV, C/(2£ 2 ) = 2100^5 (Tesla) x e 2 /(e£), 
AT 6 [350,370], and m t = 118. The Fermi energies of 
holes arc shown as the horizontal dotted lines in Fig. [S] 
For the selected parameters and magnetic field ranges, 
the bell-shape antidot is found to have the properties 
that (i) MDD ground states satisfy JVj > JVf > m t , (ii) 
only the spin-down and spin-flip transitions appear, and 
(iii) both the transitions are periodic with B (the periods 
of spin-down and spin-flip transitions are ABj = 0.0113 
T and ABk = 0.0288 T, respectively); there are no spin- 
up transition. Figure 0] displays the chemical potential 
for a stronger confinement potential, whose parameters 
are Ml = 1.9 meV, C/(2£ 2 ) = 3600 y/B (Tesla) x e 2 /(ef), 
JV 6 [368,384], and mt = 111. The differences between 
JVj and JVj are smaller than those in Fig. |SJ 

In Fig. EH we calculate the spin configuration (JVf, JVj) 
of MDD ground states and the energy difference between 
the ground and first excited states as a function of B. 
The topmost Fermi level of holes in Fig. [S] is chosen 
for these calculations. In this figure, the system starts 
in |iV T , A^.) = |128, 233) and changes into [128,234) at 
about 0.01 T larger magnetic field. This spin-down tran- 
sition |JVf, JVj) — > |JVf , JVj + 1) corresponds to the chemi- 
cal potential jumps in Fig. [S] occurring whenever each 



E = 274.16 meV 

F 




B [T] 

FIG. 9: Energy differences of the ground and first excited 
states of the antidot studied in Fig. |H| as a function of B. 
We choose the topmost Fermi level shown in Fig. [HI so that 
the transition of MDD ground states (their spin configura- 
tions are marked in this figure) follows the top zigzag line of 
Fig. |U At a cusp in Ei(B) — Eo(B), the spin-down tran- 
sition |JVf,iV|) — > |iVf,iVj + 1) occurs. The spin-flip tran- 
sitions | JV-f , iVj_ ) — * \Nf + 1, JVj. — 1) are also found around 
B = 1.21375, 1.2425, 1.27125, and 1.3T. A domain boundary 
between two MDD ground states is indicated by the numbers 
(JV|,JVj). The state right to the domain boundary has JVf 
and JVj numbers of spin- up and -down holes, respectively. 



chemical potential line intersects the Fermi level (see 
the changes in the occupation numbers of spin-up and 
-down electrons in Fig. As B increases further, there 
appears a spin-flip transition of 1 128, 234) — > 1 129, 233) 
around B = 1.215 T (see the corresponding degener- 
acy point marked by a square dot in Fig. [S] and the 
changes in the occupation numbers of spin-up and -down 
electrons in Fig. OJ. For the studied bell-shape anti- 
dots, we find only the spin-flip transition of the type 
\Nf, JVj) — > | JVf + 1, JVj — 1) with increasing B, where the 
total hole (electron) spin increases (decreases). This can 
be understood from the fact that for the bell-shape anti- 
dot with a negative potential curvature near MDD edges 
( JVj, JVj > m t ), the confinement potential behaves like 
~ B for the states with m > m t [see Eq. ©], while the 
Coulomb matrix elements goes like ~ \ B~. Consequently, 
the total confinement energy of holes increases faster than 
the total Coulomb energy as B increases. This causes the 
hole droplet to minimize its size, i.e., JVj becomes smaller 
for stronger magnetic field. This feature is absent in the 
inverse parabolic antidot, where the spin-flip transition 
is found to occur as |JV T ,JVj) -> |JV T - 1, JVj + 1) with 
increasing B (see Sec. VII). 

However, the spatial splitting (~ ~ of 

spin-up and -down edges of MDDs is much larger than I 
in the bell-shape antidot studied in Fig. [SJ In this case 
the coupling of the spin- up part of the MDDs to the ex- 
tended edge channels will be negligibly small, compared 
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to that of the spin-down part, and thus, the Kondo res- 
onance can not occur at the spin-flip transition points. 
Instead, these parameters may represent a good model of 
(h/2e) Aharonov-Bohm oscillations without Kondo reso- 
nances measured in Ref. |9j. To see the Kondo resonance, 
one should have the conditions 0] that (i) there exists 
the spin-flip transitions and (ii) both the spin-up and 
spin-down electrons are weakly coupled to extended edge 
channels; the condition (ii) is equivalent to — i?| < I 
in MDD states, where i? T ~ ^/2N-\t and R l ~ y/2N\e 
are the radii of spin-up and -down droplets, respectively. 
It is possible to make R± — i?| smaller by making the 
confinement potential stronger. Then, the corresponding 
MDD would be closer to the actual electronic state mea- 
sured in Ref. ^2] • Figure 0] displays such a case with 
smaller i?| — i£j . 

As will be shown in the next section VI, the spin-down 
and spin-flip transitions found in the bell-shape antidots 
are in good agreement with the predictions of the capac- 
itive interaction model ^5|, and thus they can explain 
qualitatively the experimental conductance data |12( of 
Aharonov-Bohm oscillations with Kondo- like signatures. 
However, for realistic quantitative comparison, the cou- 
pling to higher Landau levels by electron-electron inter- 
actions may have to be included. 



VI. CAPACITIVE INTERACTION MODEL OF 
BELL-SHAPE ANTIDOTS 



Both the spin-down and spin-flip transitions shown in 
Figs. 0]and|Blare almost periodic with B, indicating that 
the bell-shape antidot with more than 300 holes may be 
large enough to be described by the phenomenological ca- 
pacitive interaction model for antidot |l5j | . Our Hartree- 
Fock results are indeed consistent with the model. Below 
this is demonstrated by determining the parameters of 
the model from the Hartree-Fock results of the antidot 
studied in Fig. [HJ For the antidot in Fig. 0]one can get 
the parameters in the same way. 



A. Excess charges 

The capacitive interaction model was developed for the 
case that the local filling factor is two around antidot. In 
this model, the excess charges ^Ej can be defined as 



6q a (N a ,B) = e(N a -N rT (B)) 



(15) 



in terms of the number of holes N a and a =T,4- The 
function N„(B) has the meaning of "the optimal num- 
ber" of spin-er holes that minimizes the total energy in 
the absence of holes with the opposite spin. Note that 
N<j(B) is real- valued while the actual number of holes is 
an integer, which prohibits continuous change of N a and 
leads to Coulomb blockade. Of course, the detailed form 
of N a (B) depends on the shape of the antidot potential 



near its edge. For a sufficiently small range of B, one can 
use an approximately linear form 



eN a (B) =e{a a B + b a ). 



(16) 



Here, a a e is the rate of excess charge accumulation with 
increasing B and b a e originates from the positive back- 
ground charge. Both a a and b a are taken as constants. 
This is a good approximation for large-size antidots in 
strong magnetic fields. Then, one has 



Sq a {N a ,B) 



= N n - a„B - b a 



(17) 



By assuming that the excess charges interact capaci- 
tiv ely. one can write the total energy of an isolate antidot 

as m 

E C l{5q],5qi) = ^^25q a (C~ 1 ) a , a ,5q a , 



a, a' 



(5 qi +a5 qi ) 2 &q 2 } 



2C ut 2Cm 



(18) 



where C is a capacitive matrix, a = |CjjJ/Cjj, C ou t — 
Cn — a|CjJ, and Cj„ = C^. These elements of capaci- 
tive matrix can be taken as constants for large-size anti- 
dot with ABab "C B, where C CTjCT / can vary very slowly 
over several Aharonov-Bohm periods ABab- Thus, for 
large-size antidots in strong magnetic fields, the capaci- 
tive interaction model defined by Eqs. I|17I18[) is a good 
approximation. The model is determined by the constant 
parameters a, Cj n , C ou t, a a , and b a , and is analogous to 
the constant interaction model of quantum dots [25[ . 



B. Transition conditions 

The conditions of transitions between MDD ground 
states can be rewritten by using Eci- For example, the 
condition (|13fl for the spin-down transition of |iV|, Ny) — > 
l-^T'-^l ± !) becomes E C i(Sq^,5q l ± e) = Eci(%,<%), 
which is equivalent to 



e 2 



S gi (N h B) 



0. 



(19) 



It is useful to note that N± increases (decreases) with 
increasing (decreasing) B in the spin-down transitions, 
as discussed in Sec. VA. This observation makes it easy 
to choose signs in Eq. (|19fl . 

The spin-flip transition in Eq. I|14(l can be written as 
-E , ci(<5<2 , T = ' =e ' ^91 "Fe) = Ed(Sqi ,Sq^). When the transition 
|iV|,7V^) — * | TV-p ± 1 , iV_j_=F 1) appears at B, one can rewrites 
this condition as 



%A (sq l (N LT LB)+a5q 1 (N 1 ±l, B) 



-J-«5g T (iV T ±i B) = 0. (20) 

O in ^ 
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In the spin flip transition, N± decreases (increases) with 
increasing (decreasing) B. Again this observation makes 
it easy to choose signs in Eq. Il2l)[l . When a = 1, the 
spin-flip condition of Eq. 12Ufl is reduced to a simple 
form 

M = ± I. (21) 
e 2 K ' 

C. Evolution of excess charges 

As B varies, excess charges are accumulated following 
Eq. (fTTjl . The pair of excess charges (5q^ /e,8qje) thus 
evolves with satisfying 

% = ^T£i +/3j (22) 
e e a-f 

where (3 = + (&| — A^)(aj/af) — 6;. This con- 
dition describes a line in the space of {Sq^/e, Sqi/e). 
If the magnetic field B is tuned such that one of the 
spin-up, spin-down, and spin-flip transition conditions 
in Eqs. |I1 Elm is satisfied, (N h N t ) and (N\,N[) 
become degenerate. At this degenerate point the ex- 
cess charges 5q a (N ai B) jump to their transition values 
5q a (Nl 7 , B). As a result, the evolution of the excess- 
charge pair (Sq^ /e, Sq^ / e) is restricted within a hexagonal 
cell whose boundaries are determined by the six transi- 
tion conditions in Eqs. H21 C31 H3I> d5j- Note that the 
transition conditions can be easily rewritten in terms of 
5q a as in Eqs. itl^l l2T)|l . The evolution trajectory of the 
excess-charge pair in the cell can be simply drawn when 
the value of the pair at starting B is known: The trajec- 
tory evolves along a line parallel to Sqi — (ai/a^)5qt and 
jumps to the corresponding opposite boundary when it 
collides with a cell boundary. 

D. Parameters of capacitive interaction model 

One can obtain the parameters such as a, a CT , and 
b a of the capacitive interaction model from the Hartree- 
Fock results of the ^-dependence of MDD ground states 
|iV| , Ni). Below, we show such a procedure for a ground- 
state transition on the chemical potential line which has 
the value /x = 273.33 meV in Fig. [SJ 

For this case, one can find the evolution sequence of 
\Nf,Ni) as a function of B in Fig. |S] (see the fourth 
zigzag line from the top one as well as the caption). From 
the periodic spin-down and spin-flip transitions one can 
first determine a and a a . Applying Eq. (|19fl to three 
consecutive spin-down transitions in B 6 (1.205, 1.245) 
T, we get a — 1 and a-f = 88.89. Similarly, applying 
Eq. (|21|) to two consecutive spin-flip transitions at B — 
1.2125 T and 1.24125 T, we find a T = 34.78 and oj. = 
54.11. For a = 1, it is noteworthy that the spin- flip 
condition in Eq. (|20|) can be simply reduced to Eq. (|21|l , 
and that the evolution trajectory of (Sq-\ , Sqi) is restricted 



1 - 




-0.5 0.5 

5q<f/e 



FIG. 10: The evolution trajectory (solid line) of 
(Sq^(B), 8qi(B)) for the bell-shape antidot potential, while 
the MDD ground state evolves following the fourth zigzag 
solid line in Fig. |5] The conditions of the spin-down [Eq. 
l(T3jl ] and spin-flip [Eq. I|140 ] transitions define the cell bound- 
aries (dashed and dotted lines, respectively). The trajectory 
is determined by Eq. I|16jl . thus it evolves parallel to the line 
of Sqi — (ai/ai)Sqi . When the trajectory collides with a 
cell boundary, it jumps following vertical (spin-down) or di- 
agonal dashed-dot lines (spin- flip transition). The trajectory 
starts at B = 1.205 T and ends at 1.245 T, and the arrows 
represent the direction of evolution with increasing B. The 
parameters (a = 1, a;/a| = 1.55, b] — 85.33, fo| = 165.79) 
are determined from the result in Fig. |S](see text). 

in a parallelogram, instead of a hexagonal cell, since no 
spin-up transition appears. 

The other parameters b a 's can be uniquely determined 
from the assumption that the evolution speed (the accu- 
mulation speed of excess charges) is constant with vary- 
ing B and from the facts that (i) the trajectory moves 
within the cell and (ii) at the spin-flip transitions (e.g., 
at B = 1.2125T in this case) the value of Sq^ is fixed by 
Eq. (|21|l when a = 1. Combining these with Eq. i|17|) 
we obtain b^ = 85.33 and b± = 165.79. Note that one 
can not determine C ou t and C'i n from the data in Fig. |H| 
when a = 1. 



E. Evolution trajectory of excess charge 

We draw the evolution trajectory (<5gj , Sqi) by using 
Eq. (|22l) and the parameters of a, a a , and b a obtained 
in the last subsection VID. The evolution trajectory de- 
scribes well the features of the ground state transitions 
in the range B € (1.205, 1.245) T, such as the sequence 
of transitions and the differences of magnetic fields be- 
tween adjacent transitions. In the range of B > 1.25 



T, the trajectory becomes to deviate from the Hartree- 
Fock results in Fig. [8] indicating that the parameters 
are not constant but slowly vary with B. This finding 
demonstrates that the capacitive interaction model can 
describe the behavior of the bell-shape antidot very well. 

Some predictions of the capacitive interaction model 
are useful in understanding of microscopic numerical re- 
sults. It predicts that a = 1 when the interaction be- 
tween Sqi and 6q± is maximal, indicating that the spin-up 
and -down excess charges of MDD states formed in the 
bell-shape antidot are strongly coupled. It also predicts 
that when a = 1, no spin- up transition appears due to 
Coulomb blockade while periodic spin-down and spin-flip 
transitions with B can occur (the periods are AB^ and 
ABk, respectively). When a = 1, the ratio of the tran- 
sition periods relates to the ratio of accumulation speeds 
of excess charges as ABk / AB± = 1 + aja^. 

The ratio of accumulation speed a[/a-\ of excess 
charges generally relates to the ratio of hole occupation 
area R 2 /R 2 , since larger area gives smaller Aharonov- 
Bohm period and thus faster accumulation of excess 
charge. Our system has the area ratio R 2 /R 2 = N^/N^ ~ 
1.8 when (N^N^ = (128,233), while the speed ratio 
is ai/a^ = 1.56. The discrepancy between aj./a| and 
R 2 /R 2 originates from the fact that the accumulation 
speed depends on antidot potential shape. The ground 
states of the tested antidot have large spatial splitting be- 
tween the edges of spin-up and spin-down states so that 
the outmost spin-up orbital at m — N-\ sees a different 
potential slope from the spin-down orbital at m = N[. 
For the case of iVj ~ iVj we expect a[/a-\ — R 2 /R 2 . 

When aja^ ~ 1 (i.e., Aj ~ AjJ, one gets ABk — 
2AB± from the relation of AB K /ABi = l+a x /a T . Then, 
one has two spin-down and one spin-flip transitions in one 
Aharonov-Bohm period, which is the type (i) evolution 
trajectory of excess charges predicted in Ref. 15]. The 
type (i) trajectory matches very well with the experi- 
mental findings of (h/2e) AB oscillation and Kondo-likc 
signatures in Ref. |l2(. 



VII. HARTREE-FOCK RESULTS OF INVERSE 
PARABOLIC ANTIDOTS 

In this section, we consider the inverse parabolic an- 
tidot shown in Eqs. (121 If) . From the Hartree-Fock 
calculation, the ground state is found to be a MDD when 
we choose the parameters M7 = 1.5 meV and B ~ 1.45 
T. The transition of MDD ground states in this anti- 
dot is plotted in Fig. ^] Three types of transitions 
can appear. As B increases, while the direction of the 
spin-up (-down) transition is |A|, Aj) — > |jVf + 1,N±) 
(|JVf , JVj_) — > l-ZV-f , Ni + 1)), the spin-flip transition occurs 

as |A T , A x ) l-^T ~ !> N l + !)• Note that the direc- 
tion of the spin-flip transition is opposite to the bell- 
shape case. This can be understood from the fact that 
~ 1/B in the inverse parabolic case [see Eq. (JSJ], while 
Coulomb matrix elements goes like ~ \[B. Due to this 
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FIG. 11: Magnetic field dependence of the chemical poten- 
tial for an inverse parabolic antidot with Ml = 1.5 meV. 
Chemical potentials are shown for two different values of hole 
Fermi energies. Squares indicate spin-flip transitions, while 
vertical jumps show spin-up or spin-down transitions. In 
the spin-flip transitions, the occupation number configura- 
tion changes from the lower to the upper panels of Fig. [3] 
as B increases. For the Fermi level drawn by the dashed 
horizontal line, | iV-f , iVj_ ) evolves as 1 66, 67) =>• |65, 68) — > 
|66, 68} |65, 69) |64, 70) |63,71) |62, 72) as B 
increases from 1.4 T. Here =>■ indicates the spin-flip transi- 
tion, while — » shows the spin-up or -down transition. For the 
Fermi level drawn by the solid horizontal line, |7Vj , jVj,) evolves 
as 66, 67) => |65, 68) -» |65, 69) =>■ |64, 70) => |63, 71) 
1 62 , 72). The main different result between the upper and 
lower Fermi energies is that the vertical jumps represent the 
spin-up and -down transitions, respectively. 



magnetic field dependence, the arguments based on the 
total energy minimization indicates that the hole droplet 
can have a larger size (larger AjJ for a stronger B. The 
features of spin-flip transition around an antidot thus de- 
pend on its potential shape. The Kondo effect associated 
with similar spin-flip transitions in a parabolic electron 
quantum dot is studied theoretically in Ref. [2(|. 



In the case of the inverse parabolic antidot, the spin- 
flip transitions cause that N± becomes larger with B. 
And, they appear more frequently with varying B than 
the spin- up and spin-down transitions. As a result, the 
hole spin polarization (Nf — Ni) decreases as B increases. 
In electron language, this means that the electron spin 
polarization around the antidot increases with B. These 
features cannot be explained by the predictions of the ca- 
pacitive interaction model 0] and do not match with the 
experimental data of Kataoka et al which indicate 
that the spin polarization around the antidot is periodic 
with B. It would be interesting to investigate the prop- 
erties of an inverse parabolic antidot experimentally. 
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VIII. CONCLUDING REMARKS 

We have investigated electronic properties of antidots 
in the integer quantum Hall regime (the local filling fac- 
tor around the antidot is two), by using an electron-hole 
transformation and a Hartree-Fock approach. Our nu- 
merical work shows that when the antidot potential is 
strong enough, hole MDDs of a single-Slater-determinant 
form are stable in a certain parameter range and repre- 
sent exactly the ground states. 

For a bell-shape antidot with more than 300 holes, we 
find that there exist the spin-down and spin-flip tran- 
sitions between MDD ground states as magnetic field 
varies, and that their properties, such as sequence of tran- 
sitions and the spin polarization of MDDs, are in good 
agreement with the phcnomenological capacitive interac- 
tion model. Thus, hole MDDs and the bell-shape antidot 
could be a good model for the antidot studied experimen- 
tally by Kataoka et al. The Hartree-Fock results of the 
ground-state transitions allow us to obtain the parame- 
ters of the capacitive interaction model and the accumu- 
lation of excess charges around the antidot as a function 
of magnetic field. 

The properties of hole MDDs depend on the competi- 
tion between electron-electron interactions and hole con- 
finement potential. The competition determines the de- 
gree of spin splitting between spin-up and -down edges 
of MDD: stronger hole confinement (weaker Coulomb en- 
ergy) favors the droplets with smaller size. As a conse- 
quence of this competition the direction of the spin-flip 



transitions depends on the detailed potential shape (e.g., 
the potential curvature near MDD edges). For example, 
in the bell-shape antidot the spin-flip transition occurs 
as |iV T , Ni) -> \N-[ + 1, N± - 1) with increasing B, while 
|JVf, Ni) — » |iVf — 1, N± + 1) in the inverse parabolic an- 
tidot. Thus, the physical properties of antidots are con- 
nected to the shape of the antidot potential. It would be 
interesting to investigate antidots with different shapes 
and compare their properties with each other. 

For a better quantitative comparison with the experi- 
mental data the coupling to higher Landau levels by the 
Coulomb scattering may have to be included. However, 
it is a challenge to calculate the matrix elements of the 
screened Coulomb interaction analytically. It is also de- 
sirable to use a more smoothly varying antidot potential 
than the one we used here (the first derivative our po- 
tential as a function of m changes suddenly at mt). 

Finally, reconstruction of the edge structure around 
the antidot may take place, similar to the extended edge 
channels [2?], I2§1 129| . The electron-hole transformation 
makes it possible to investigate numerically the edge 
reconstruction near the antidot, which requires a proper 
treatment of correlation effects beyond Hartree-Fock 
approximation^^ [3(j ■ 
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